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ABSTRACT 

We develop a simple model for computing planetary formation based on the core 
instability model for the gas accretion and the oligarchic growth regime for the ac- 
cretion of the solid core. In this model several planets can form simultaneously in the 
disc, a fact that has important implications specially for the changes in the dynamic of 
the planetesimals and the growth of the cores since we consider the collision between 
them as a source of potential growth. The type I and II migration of the embryos 
and the migration of the planetesimals due to the interaction with the disc of gas 
are also taken into account. With this model we consider different initial conditions 
to generate a variety of planetary systems and analyse them statistically. We explore 
the effects of using different type I migration rates on the final number of planets 
formed per planetary system such as on the distribution of masses and semimajor 
axis of extrasolar planets, where we also analyse the implications of considering differ- 
ent gas accretion rates. A particularly interesting result is the generation of a larger 
population of habitable planets when the gas accretion rate and type I migration are 
slower. 
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1 INTRODUCTION 

Planetary astronomy is a young science, and until recently, 
1 was essentially devoted to the study of planetary bodies in 
our own Solar System. The extrasolar planets found since 
1995 have vastly expanded our database by increasing the 
number of known planets by more than 200. Although the 
distribution of masses and semi-major axis of observed ex- 
trasolar planets is highly biased towards those planets that 
are detectable using Doppler radial velocity technique, the 
increase in precision and continuity of the surveys has given 
the possibility of unveiling a large variety of planets and re- 
formulate the theories of planetary formation. As the num- 
ber of detections grows, statistical studies of the properties 
of exoplanets and their host stars can be conducted to un- 
ravel some of the processes leading to the formation of plan- 
etary systems. 

In this frame several models o f planetary form ation 
have been presented in the last years. IChambera l|2006l ) and 
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iThommes et al.l (|2003l ) have developed semi analytical mod- 
els based on the oligarchic growth regime, that retain much 
of the simplicity of the analytic mo dels, but incorpora te 
more physics from the detailed ones. Ilda fc LinI (|2004al lbl. 
l2008h have studied the mass and semimajor axis distribu- 
tion of extrasolar planets through a very simple model but 
considering all the relevant physics involved in the process 
of planetary formation. 

We begin our s tudies in a previous paper 
l|Miguel fc Bruninil |2008| ). where we developed a very 
simple model based on the core instability model. The 
accretion of solids was based in the regime where the 
largest embryos dominate the dynami cs of the planetesimal 
disc: the oligarchic growth regime (|Kokubo fc Idal Il998l : 
Ilda fc Makinolll993l ). When the cores are large enough to 
have an associate envelope which can not be sustained by 
the hydrostatic equilibrium any more, the gas accretion 
process begins. In our previous work we explored the effects 
of different gas accretion rates on the final mass distribution 
of extrasolar planets, and found a strong dependence. But 
our simple model allowed us to form only one core per disc 
and did not allow us to consider the dynamical evolution of 
the cores or the planetesimal disc. 

In this paper, we have improved the model introduced 
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in Paper I considering other important effects. The mainly 
improvements are the formation and evolution of several 
cores simultaneously in the protoplanetary nebula, a fact 
that also influences the dynamic of the planetesimal disc. 
The embryos can grow by the accretion of planetesimals or 
by accreting other embryos, in this work we consider the 
collision among them as a source of growth. 

We also incorporate the orbital evolution of the pro- 
toplanets due to its tidal interaction with the gas in the 
surrounding protoplanetary nebula, considering two regimes 
of planetary migration: type I and II regimes, which result 
mostly in a radial motion towards the host star. Type I 
migration involves low mass protoplanets who are rapidly 
moved towards the star l|Goldreich fc Tremaine|[l980l : IWardI 
1 19861 ) and type II is the regime of migration for those pro- 
toplanets whose mass is large enough to allow them to open 
up a gap in its orbit and migrate on a rate contr olled by the 
disc viscous time-scale (|Lin fc Papaloizoulll985^ . This pro- 
cesses could have important consequences on the final mass 
and semimajor axis distribution, such as on the final number 
of planets formed per protoplanetary disc. 

The disc of planetesimals also interacts with the nebu- 
lar gas. This gas drag effect causes the migration of plan- 
etesimals towards the host star before they become large 
enough to decouple from the disc gas. This effect, which 
was included in our model, affects directly the distribution 
of solids in the disc and the growth of the cores. 

We show that a variety of planetary systems can be 
constructed using these principles, varying the initial condi- 
tions, and that the mass and semimajor distribution of ex- 
trasolar planets is strongly dependent on the gas accretion 
model and type I migration rate considered, being larger the 
population of habitable planets when the rate of gas accre- 
tion and type I migration are slower. 



2 DESCRIPTION OF THE MODEL 

In our previous work (jMiguel fc Bruni"nill2008f ). Paper I, we 
developed a simple mo del for computing planetary forma- 
tion based in a work of llda fc LinI (|2004al ) (hereafter IL04), 
but in corporating the accretion rates given bv lFortier et al.. 
l|2007l ). This lead us to a very different results regarding 
the final mass distribution of extrasolar planets, which was 
found strongly dependent on the gas accretion model con- 
sidered. 

With the aim of improve the model introduced in Paper 
I, we have extended it by considering other important effects. 
In this section we explain briefiy our previous model and 
show the improvements made. 

2.1 The previous model 

For the sake of completeness, we will explain our previous 
model. Although very simple, it incorporates all the impor- 
tant aspects of planetary formation. 

We considered one initial core per protoplanetary disc, 
with an initial mass of lO^^M®. The in itial disc was b ased 
on the minimum mass solar nebula of iHavashil (|l98ll )rfor 
more details see section 12. 2p . and it was defined between 
0.01 and 100a«. The population of the nebular mass scal- 
ing parameter, fd, which states the solid mass in the disc 



in terms of the minimum mass solar nebula model, was dis- 
tributed with a Gaussian distribution in terms of logw{fd) 
with dispersion 1 and centred at 0.25, in order to charac- 
terise different discs. We supposed that fg, which is a scaling 
parameter for the gas disc, was equal to fd, this correspond 
to discs with solar metallicity. 

The solid and gaseous disc were not time-invariant. The 
gaseous disc changed globally, decaying exponentially with 
a characteristic time-scale of Axlif years, in accordance to 
curre nt estimates of disc lifetim es around young solar type 
stars (|Beckwith fc Sargent|[l99^ and the solid disc changed 
locally, suffering the depletion of planetesimals produced by 
the effect of core's accretion. 

The lO'^^M© initial core was located with equal prob- 
ability per interval of log{a), and it grew with the fol- 
lowing solids accretion rate based on the particle-in-a-box- 
approximation (Safronov 1969), 

where R is the physical radius of the core, Qk is the Kepler 
frequency, Mt is the total mass of the core, Ed is the solids 
disc surface density and a is the relative velocity between 
the embryo and the disc of planetesimals. 

The factor 10.5 3 was introduced for consideri ng the F 
factor introduced bv lGreenzweig fc Lissauerl (|l992l ). and the 
approximations made for the eccentricity e, the inclination 
i and the disc scale of high h{a) in the high-cr equilibrium 
regime, 2i ~ e and h{a) ~ ai. 

Finally, a factor of 4 was introduced in the above equa- 
tion, which is an approximate value taken in orde r to fit 
the solid accretion rate used bv lFortier et al.l (|2007l ). which 
includes the evolution of the planetesimals r.m.s. e and i 
and the drag effect caused by the gaseous envelope on the 
incoming planetesimals. When the gravitational perturba- 
tion due to the protoplanets is balanced by dissipation due 
to the gas drag, the planetesimals r.m.s. ecc entricity attains 
an equilibrium value, which was obtained bv lThommes et all 
(2003). Th is effect was c onside red in the simulations per- 
formed by iFortier et all (|2007l ). where the solid accretion 
rate was prescribed as that obtained by ommes et al.l 
(|2003l ). and where the enhancement of the protoplanet's 
cross section due to the drag effect caused by the proto- 
planet's gaseous envelope on the entering planetesimals was 
also considered. We compared the core accretion rates found 
without considering these eff ects with all the cases consid- 
ered bv iFortier et al.l ()2007l ). The results showed that the 
core accretion rate is ^ 4 ti mes larger when the evo lution of 
the e and i are considered (jMiguel fc Bruninlll2008l ). 

When a core of mass Mc reached a certain critical mass, 
it started the gas accretion on to the core. The critical 
mass at which the gas accretion process began was obtained 
analytically by ^St evenso n (1982) and it was improved by 
llkoma et al.l l|200ol ) through numerical simulations. We con- 
sidered a simplified version, 

Mc„t ~ (ttt^ t) ' (2) 

The gas accretion rate was obtained by fitting the re- 
sults of the self-consistent code developed bv IFortier et al.l 
|2003) 
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dt 



Mt 



where Mg is the mass of the surrounding envelope and Tg is 
its characteristic growth time, 



1. 64x10 



9( Mt 



yrs 



(4) 



This process stopped when the protoplanet consumed 
all the gas available on its feeding zone, or what is the same, 
when it reached its isolated gas mass, given by. 



MgAao = 27ra6ri/E„ 



(5) 



where Eg is the gas surface density and brn is the typical 
distance betwe e n two adjacent embryos, which according to 
iKokubo fc Jdi I 1998l l is ~ WrH, being rn the Hill radius. 



2.2 The protoplanetary nebula and the initial 
cores 



As in our previous paper l|Miguel fc Bruninil |2008| ). the 
structure of the protoplanetary nebu la is ba sed on the min- 
imum mass solar nebula (MMSN) of Hayashi (1981j), where 
the surface density of solids is defined as follows: 



a 
lau 



gem 



(6) 



where rjice is 1 inside the ice condensation radius and 4 out- 
side it, expressing the effect of water ice formation across the 
ice boundary. The "snow line" is located at aice ~ 2.7(-^|-) 
au from the central star of mass M*. 

Similarly, the volume d ensity of gas fol l ows a n expo- 
nential profile, as explained iThommes et al.l (|2003h . which 
is 



p(a, z) = pg,o(a)c''("'^ gem 

where pg^ is the mid-plane value density given by. 



Pg,o{a) = 1.43:10-^9 



a 
lau 



gem 



(7) 



(8) 



Unlike the previous calculations, we assume that the 
surface density of gas everywhere varies with time as: 



Pg = Pg.oe 



(9) 



with Tdisc the gas depletion time-scale. 

We consider a disc between Rin and 30 au, where the 
inner disc radius is estimated by. 



R,„ = 0.0344* 



1500° iC 

Tsuh 



(10) 



as explained IVinkovid l|2006l ). In the above equation '3/ is 
a factor of ~ 2 which depends on the disc structure and 
radiative transfer, the dust sublimation temperature, Tgub, 
is taken as 1500° A" and L* and Lq are the stellar and Sun 
luminosity respectively. 

In this work we consider stars with masses between 0.7 
and 1.4M0, the stellar luminosity to this range of masses is 
given by 



(3) I^^ Ml. 



Me 



(11) 



In our previous work we considered one core per disc, 
here we form several cores in the same disc. Initially we start 
the simulation with a number Npia„ets.o of cores through the 
disc, separated by lOrn- This could have important conse- 
quences on the final distribution of masses and semimajor 
axis of extrasolar planets, specially for the changes in the 
dynamic of the planetesimals, and also considering the colli- 
sions between em bryos as a very important sou rce of growth, 
as was shown by iBrunini fc Benvenutj l|2008t ). 

In the core instability model the first step in the plane- 
tary formation is the coagulation of dust particles. When 
particles of many meters in size have been formed, the 
gravitational forces became more important than gas drag, 
and the larger planetesimals grow faster than the smaller 
ones. This st age of increasingly rapid gro wth is the run - 
away growth (Creenberg et aLlll978 : Kokubo fc Idalll996l ). 
In this stage, the planetesimal dynamic is dominated by the 
interaction with other planetesimals, until the larger ones 
grows so much and start t o dominate the veloc ity distribu- 
tion of the planetesimals, l|Kokubo fc Idairi998l ). This is the 
oligarchic growth stage, and is the stage where the embryos 
form. In this high-a regime, collisions between planetesimals 
tend to be destructive leaving fragments of planetesimals as 
a result of the impacts. This fragments are small enough to 
be rapidly damped by the gas drag and as a consequence, 
gravitational focusing favours the embryo, increasing its ac- 
cretion rate (|Chambersll2006l ). This could be an important 
factor in the growth of the embryos, but is not taken into 
acc ount in this wor k. 

llda fc Makiiiol (Il993l ) found through N-body simula- 
tions, the minimum mass necessary for starting the oli- 
garchic growth stage, which is given by 



Mali 



(12) 



where m is the effective planetesimal mass, which was found 
adopting the next a cumulative power law mass distribution 
of planetesimals. 



N{m)dm oc m ^dm 



(13) 



where p is 1.8 — 3. We consider p — 2 . 5, consistent with the 
resuhs of iKokubo fc Idal (|l998l , ]2000l ) for the spectrum of 
masses of a population that has relaxed to isolated runaway 
bodies. Then the typical mass we use is found through 



m'N{m')dm' 



(14) 



here we suppose mmin ~ 6.Sxl0^^g and mmax = 6.3xl0'^^g 
are the minimum and maximum planetesimal mass, equiv- 
ale nt to 0.1 — lOOKm ra,dius (e .g. the same values adopted 
bv lBrunini fc Benvenutd l|2008l ')'). here we note that the disc 
mass is contained in the small bodies. The typical planetes- 
imal mass found with equation (|14p is the one used and is a 
constant in all the simulations performed. 

One initial core is located at a = Rin, the rest of the 
cores are separated 10 rn each other until the end of the 
disc is reached. Their initial masses are given by equation 
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l|12p . where we see that different location in the disc leads 
to a different initial oligarchic mass. 

2.3 The growth of the cores 

The solid accretion rate for a core in the oligarchic growth 
regime, conside ring the particle-in-a-box approximation 



Safronovl [19691 ') is given by equation ([T]). iThommes et alj 



20031 ) found a expression considering the evolution of the 



planetesimal rms e and i. They found that when the em- 
bryos' gravitational perturbation are balanced by dissipa- 
tion due to the gas drag, the planetesimal rms e attains an 
equilibrium value given by. 



1.7 ■ TUTS 



(15) 



where Cd is a dimensionless drag coefficient which is of order 
unity, and pm is the planetesimal bulk density, which is ~ 
l.Sgcm^'^. The equation ([l]) could be rewritten using this 
expression. 



dt 



3.96^C^G5Af,V|„, 



Pm PliO. If m 15 



-Ml 



(16) 



where pM is the embryo bulk density (pjv/ ~ 1.5(;cm~^). 

The main difference between equation (|16p and the one 
considered in Paper I is that in the above equation we take 
into account the evolution of rms e and i through the equa- 
tion (lisp , and in Paper I we ap proximated this evo lution 
considering the cases analysed bv lFortier et al.l (2007). 

Equilibrium rms velocities are a good approximation as 
long as the embryos reach roughly one Earth m ass, a state 
that is quickly attained in the regions a < lOau (jChambersI 
I2OO6I ). Nevertheless, our model intent to remain as simple 
as possible, in such a way that we can perform statistics of 
a large number of simulated result. Within this in mind, a 
more appropriate velocity evolution model is for the moment 
out of our possibilities. 

The growth of the cores terminate when they consume 
all the planetesimals on their feeding zone (Aac = brn), 
or equivalently when the solid superficial density is Ed = 
(see section [2^ . But they also can stop the growth when the 
density of planetesimals is substantially depleted by ejection 
((rhommcs ct al. 2003; Ida & Lin 20043), this is, when the 
ratio of collision to ejection probabilities is 



< 1 



(17) 



where — ^ 2GMi, j a is the escape velocity from the pri- 
mary and Vs — \J GMt / R is their characteristic surface ve- 
locity. 

Once the core became massive enough to retain a gas 
envelope, the effect of this atmospheric gas drag on the plan- 
etesimals increases the collision cross section of the proto- 
planet. Considering the model for a purely rad ia tive atmo- 
sphere IStevensonl Il982l : llnaba fc IkomaH2003l ). IChambersI 
l|2006t ) found an approximate expression for the enhanced 
collision radius {Re) of the embryo. 



' Rc\ _ 0.000344^^ eP / Mt 
R I ^ Kr^Ed I 



24ei 



(18) 



where c is the velocity of light, P is the orbital period of 
the protoplanet, k is the opacity of the atmosphere which is 
considered as ~ 4cm^g~^, rm is the planetesimal's typical 
radius and the equilibrium eccentricity is considered as « 2 
in this expression. 

Once the core reaches the critical mass given by equa- 
tion ©J the gas accretion process begins. The model for the 
gas accretion on to the core is the same used in Paper I, and 
it was already explain in section [2Tl1 

We also perform some simulations considering the gas 
accretion rate given by IL04, 



Tg,(IL04.) 



= l.xW 



Mt 

Mr. 



yrs 



(19) 



which differs mainly in the exponent, taken as 3 by 
them. It leads as to a smaller gas accretion rate, which has 
an enormous influence on the distribution of masses, as was 
explained on Paper I. 

Until this point, we have assumed that embryos grow by 
accreting gas and planetesimals only. However, embryos can 
also grow by accreting other embryos. When two protoplan- 
ets are too close to each other, their mutual gravitational 
perturbation induce high eccentricities, which enable their 
orbits to cross. This process may lead to close gravitational 
encounters and violent collisions between the embryos. 

We suppose that collision between protoplanets will oc- 
cur if their orbital spacing is less than 3.5 Hill radius. Con- 
sidering that the evolution of planetary atmosphere after a 
merger is a complex and poorly understood process, we anal- 
yse two different and extreme scenarios. On the one hand we 
assume that all embryo-embryo collisions will lead to coales- 
cence to form a single body where the result is simply the 
sum of the masses (gas and solids) of both embryos, no mat- 
ter the composition of the protoplanets. On the other hand 
impacts so energetic can cause a net loss of atmosphere, be- 
cause the kinetic energy of the impact is partly transferred 
to the associate envelope accelerating its molecules to veloc- 
ities greater than the escape velocity from the protoplanet. 
So we want to explore this possibility and in the other sce- 
nario, if both colliding planets have an associate envelope, 
after the impact their atmospheres would be lost and the 
result would be a merger of their cores. The mass and semi- 
major axis distribution were similar in both cases, for this 
reason we stayed with scenario 1) to not add greater com- 
plexity to the model. 



2.4 Evolution of the solids surface density 

Protoplanet ary discs are not quite steady, so the surface 
density of solids is not a time- invariant. It will change over 
time due to different effects, one of them is the depletion 
of planetesimals produced by the effect of cores' accretion 
(this effect was also considered in Paper I). If, T.d,o is the 
initial surface density, we have to remove what already ate 
the embryo, so the evolution of Ed is given by. 



Ed — Ed 



M, 



2iiahrH 



(20) 



On the other hand, the disc of planetesimals also inter- 
acts with the nebular gas. This gas drag effect will cause 
a radial motion of planetesimals before they become large 
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enough to deco uple from the disc gas. The orbital deca y 
occur on a rate lAdachi et al.lll976l : Frhommes et al.ll2003l ). 



da 



^ 49 2 1 .2 

''+16^ +8* 



(21) 



where Tgas is the characteristic time-scale of the gas drag 
and ri is the fractional deviation of the gas velocity Vgas 
from Keplerian velocity vk due to the pressure gradient, 

2 



V -- 



Vk 



VK 



IStt 
"64~ 



\VK / 



(22) 



with Cs being the sound speed and Cs/vk — h{a)/a. 

The characteristic time-scale of the gas drag is given by, 

^'""^ " {CD/2)nrlp,anK ^^^^ 

Finally, applying continuity to the planetesimal disc, 
the gas drag effect acting on the planetesimals leads to a 
rate of change of the surface density given by. 



_ _ld_f^ da' 

dt a da\ dt / 



(24) 



Equation (|20|) is added to the equation (|24|) as a sink 
term, so we obtain the total change of the solid surface den- 
sity in the disc. 



with q = 1 — 1.5, we use 1.5. 

The time-scale of this inward migration is smaller 
than the disc lifetime, since it seems very unlikely, re- 
cent works dealing with planet-disc interactions have pro- 
posed mechanisms that could slow down or stop type I 
migration. Several mechanisms have been proposed like 
the s tudy of the magneto-hydro dynamic turbulence in the 
disc (iNelson fc Papaloizoull2004l ). the effects of a disc con- 
taining a toroidal magnetic field (lTeraueml | 2003l). the tri- 
dimensional effects studied by iMasset et al. ( 20061 ) , varia- 
tion in th e tempera ture gradient and surface density in the 
disc (Ma sset et al. |[2006.), the effect of including a proper 
energy balance on the interaction of a low-mass planet with 
a protoplanetary disc (jPaardekooper fc Mellernall2006l ) . and 
so on. 

Halting or slowing down the inward migration requires 
careful computation of tidal effects between the core and the 
star, this is much beyond the capabilities of our model. For 
this reason a factor of ■ 



is introduced in equation (|25 



Cmigl 

to consider non linear effects without introduce a mayor de- 
gree of complexity to our model. If we want to slow down 
migration rates, Cmigi must be ^ 1. We perform simula- 
tions with Cmigi = 1,0.1,0.01. We also consider that the 
migration mechanism stops when the core reach the inner 
edge of the disc given by equation (|10p . 



2.5 Model for planetary migration 

An important contribution to our understanding about plan- 
etary formation and evolution was the discovery of the firsts 
Jupiter-like planets orbiting very close to its central star 
(hot-Jupiters). The difficulties associated with the forma- 
tion of this objects are most pronounced as we are closer to 
their host star, and they have awaken an interest in theories 
for the migration of protoplanetary embryos due to grav- 
itational interaction with the disc. In this section we will 
consider two regimes of planetary migration: the type I and 
II regimes, which result mostly in a radial motion towards 
the star and concern planets of low and high mass respec- 
tively. 

2. 5. 1 type I migration 

Type I migration acts on low-mass protoplanets, that is, for 
those embryos who are unable to open a gap on its orbit. 
Our model fo r this regime of planetary migration is simi- 
lar to that of llda fc LinI l|2008f ). who describe in detail the 
assumptions made. 

This process is caused by an imbalance in the tidal 
torques from inner and outer disc leading to angular mo- 
mentum exchange and making the planet drifts relative to 
the disc material. The time-scales involved in this process 
are smaller th an the discs lifetimes, and were calculated by 
iTanaka et al. l (200^) through 3D linear simulations, 

a 

a 

2.7 + l.lgVaf7A'/ Mp\2hpga'^ ) ^ 
- 2.7+l.lg/,U/e^ ^laJ \M^) ^^^^ 



2.5.2 type n migration 

According to the observed range of orbital semimajor axes, 
many giant planets are found to a very close proximity to 
their parent stars. At this distance the equilibrium tempera- 
ture is above the 170° -ft', the temperature at which ice (and 
the icy cores of giant planets) exists. Actually, giant plan- 
ets would form, preferably, at or outside of the ice line, this 
implies that they could have sufi^ered great orbital changes 
after formation. 

The model of giant planet type II migration presented 
here is a v e ry simple on e, and is essentially the same use by 
llda fc LinI (|2004al . [20o3) . so we will describe it briefiy. 

The process start s when the giant p l anet f orm a gap in 
the disc. According to lLin fc Papaloizoul (| 19931 ). a necessary 
condition to open a gap is 

^ /i (26) 

The above equation leads to a mass condition, when the 
mass of the protoplanet is larger than 

Mgar. = 3.75 ■ 10"* • f-^) ' (27) 

a gap is opened in the disc. When this happens, there are 
tidal forces acting between the protoplanet and the inner 
and outer edges of the gap, whose balance (or imbalance) 
depends on their difference in density. The imb alance can 
result on a migration of the protoplanet 's orbit (|Lin et al.l 
Il996l ) inward or outward. The migration time-scale is 

...,„^o.8.ioV.-'(^|)(t)(l^)lT^)*»"<») 

where a is a dimensionless parameter which charac- 
terises the viscosity and is taken as 10~^ through this work. 

The protoplanet will migrate towards the star if a < Rm 
and away if a > Rm, where Rm is the radius of maximum 
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viscous couple. This radius depends on the distribution of 
gas surface density, so it changes because Eg is not a con- 
stant. From the conservation of angular momentum, and 
considering the total disc mass as a constant, it can be found 
that 



Rn 



lOe' 



(29) 



Various stopping mechanisms have been proposed, but 
none of them seem to be effective in halting migration, for 
this reason we stop migration arbitrarily when the planet 
reach the inner radius of the disc, which is given by equation 

m- 

When a giant planet is migrating inwards, it will per- 
turbate the cores located on its path towards the host star. 
This encounter can cause the ejection or the accretion of the 
core by the giant planet, or neither, which means that the 
core will survive the giant planet path. 

The effect of the migrating giant plan ets on the for- 
matio n of terrestrial planets was studied bv lFogg fc Nelso"nl 
who present results of their N-body simula- 
tions of giant planet migration through an inner proto- 
planet/planetesimal disc, where the migration of the low- 
mass embryos (type I) was also considered. 

At the end of the simulations, they found that ~ 71% of 
the initial solid disc survive after the migration of the giant 
planet, a very small percentage is ejected into a hyperbolic 
orbit (less that 0.2%), and the rest is accreted. 

We introduce this results in our simulations, with a very 
simplified model where in a close encounter between a giant 
migrating inwards and another core, the last one has the 29% 
of possibilities of being accreted by the giant and the rest 
71% of surviving the giant passage with its orbit unchanged. 
We note that this is a very simple model, but in reality this 
fraction will change depending on the distance of embryos 
from the star. 

This effect does not significantly change the results, 
since we have more than one migrating giant planet in the 
disc. The passage of many giant planets will reduce the sur- 
viving probabilities of the cores, until they were finally ac- 
creted. 



3 RESULTS 

In this section we first analyse the characteristics of the 
planetary systems formed, and then discuss the effects of 
different prescriptions for the gas accretion rate and differ- 
ent retardation constants for type I migration on the mass 
and semimajor axis distribution of extrasolar planets. 



3.1 Characteristics of the planetary systems 
formed 

In each simulation we generate 1000 discs. Each system 
evolves for 10^ years, in a rich-metal disc, where the time- 
scale for the depletion of the gas, has a uniform log distribu- 
tion between 10^ and 10^ years, and the stellar mass has a 
uniform distribution in log scale in the range of 0.7 — 1.4M0. 
The initial number of planets per disc, Npianets,o, depends 
on the initial cores' mass (eauation ll2|l . the inner disc radius 



(equation [To]) and the mass of the host star. The final num- 
ber of planets, Npianets, is a result and shows the evolution 
of the planetary system. 

Each disc of gas is defined by fg,o, which has a log 
normal distribution with a dispersion of 1 and centred on 
. fq.o = 1, and the disc of solids is taken as — 10°'^/g,o 
dlda fc Linll2004bl ). This is slightly different from what we 
had considered before, because in Paper I we used discs with 
solar metallicity and here we suppose more metallic discs. 
The final number of planets per disc also depends on /g,o 
and /d, as seen in Fig. [1] where /d is plot against Npianets, 
considering different values for Cmigi- 
Figure 



la 



was obtained with G 



migl 



1, Fig. 1(b) 



considers C, 



migl 



0.1, the value taken in Fig. 1(c) was 0.01 



and the type I migration was not considered in Fig. 1(d) 



As seen in the Figures, planetary systems with a large 
number of planets {Npianets > 100) correspond to small val- 
ues of /d (< 3), and this fact does not change consider- 
ing different rates of type I migration. These discs have low 
mass and as a consequence the initial cores will not grow 
too much. When the mass of the cores is less than « O.IM®, 
they are not so affected by migration (which would favor the 
merging) and therefore the discs remain with many small 
planets, no matter the rate of migration considered. 

On the other hand the values of fd required to form 
planetary systems with few planets change when different 
migration rates are considered. When Cmigi ~ 1, discs with 
fd > 0.5 can form planetary systems with few planets. This 
can be explained with the rapid type I migration rate, which 
moves the biggest cores inwards accreting the other cores to 
its path, forming planetary systems with few planets near 
to the host star. Small values of fd implies discs less mas- 
sive, which will form planetary systems with few (terrestrial) 
planets. 

As migration slows down, the range of values of fd re- 
quired to obtain planetary systems with few planets becomes 



smaller, as seen in figures 1(b) 1(c) and 1(d) In the last fig- 
ure, when Cmigi = 0, we note that those planetary systems 
with few planets correspond to fd > 10, which leads to high 
mass discs which will generate giant planets who accrete all 
the other cores to its path towards the central star. 

We also analyse the kind of planets formed per disc. We 
consider that terrestrial planets are those with Mt < 7M^, 
giants with low percentage of gas are suppose to be the 
planets with Mt > 7Me and a percentage of gas mass less 
than < 15%, and the others, which have a larger percentage 
of gas are gas giants. The histogram presented in Fig. [2] 
shows the number of terrestrial planets (dotted line with 
filled circles), giants with low percentage of gas (dashed line 
with filled triangles) and gas giants (solid line with squares) 
formed per planetary system. 

Figure [2(a) | is the one obtained when Cmtgi = 1. We 
can see that planetary systems with few planets have mostly 
terrestrial ones and there are few with giant planets. As we 
explained, for this rapid migration rate, the range of discs 
which form planetary systems with low Npianets are those 
with 0.5 < fd < 30. In this range most of the discs have 
low mass and cores have less time to grow because they are 
rapidly moved to the central star. As a result most of the 
planetary systems with few planets will form only terrestrial 
planets. On the other hand those discs with a large value of 
fd are able to form planetary systems with giant planets, but 
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Figure 1. The final number of planets is plot against the disc 
mass scaling parameter, f^. Fig. |l(a)| shows the results when 
Cmigl = 1 is considered, fig. |l(b)| and |l(c)| present the results 
when Cjnigl is equal to 0.1 and 0.01 respectively and in Fig. |l(d)] 
Cmigl is which means that type I migration is not considered. 
As seen in the figures, planetary systems with a large number of 
planets are those with small values of and a large disc mass 
is necessary to form planetary systems with few planets, but the 
range depends on the type I migration rate considered. 



100 



100 




100 



100 



100 



100 



Figure 2. Histogram of the final number of terrestrial (dotted 
line with filled circles), giant with low percentage of gas (dashed 
line with filled triangles) and gas giant planets (solid line with 
squares) per planetary disc. This histograms were obtained with 



different values of C, 



migl 



In Fij 



inigl 



1 



figure [2(b)] was 



obtained when Cmigl = 0.1 and in figure [2(c) [ the value of Cmigl 
is 0.01. The number of planetary systems with few terrestrial 
planets are larger when Cmig is 1, and it decreases when the 
type I migration rate slows down. On the contrary, the number of 
planetary systems with few gas giants are smaller with Cmigl = 1; 
and it increases with small values of Cmigl- 



are the less ones. We also note a big absence of planetary 
systems with an intermediate number of planets {Npianets > 
30), this is because of the rapid migration rate, which moves 
quickly the cores to the inner disc limit, so they accrete the 
other cores to its path, forming planetary systems with few 
planets. 
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When Cmigi is smaller (Fig. 2(b) with Cmigi =0.1 and 
2(c) with Cmigi = 0.01), the population of planetary sys- 



tems with few gas giant planets increases, meaning that the 
formation of gas giants is favoured with slower type I migra- 
tion rates. On the other hand, the number of giant planets 
with low percentage of gas remains low, this population is 
not affected by the rate of migration. We also note that the 
number of planetary systems with an intermediate number 
of planets increases, so when migration slows down the fi- 
nal amount of planets per disc tends to follow a uniform 
distribution. 



3.2 Mass and semimajor axis distribution 

To investigate the mass and semimajor axis of extrasolar 
planets we perform a series of numerical simulations consid- 
ering different prescriptions for the gas accretion rate and 
different retardation constants for type I migration. 

Results are displayed in Fig. |3Ia,b,c,d) and |Ha,b,c,d). 
The first ones were found by consideri ng the gas ac c retion 
rates obtained by fitting the results of iFortier et al.l (|2007l ) 
and the second ones were obtained using the gas accretion 



rates given by IL04. In the figures 3(a) and 4(a) the factor 
considered for delaying the type I migration rate is Cmigi = 
1, in Fig(3(b)l and |4(Eyi C™gj = 0.1 Fig. [3(c)] and \^c) 



those obtained with C, 



migl 



0.01 and finally Fig. 3(d) and 



4(d) where found without considering the effects of type I 
migration {Cmigi = 0). 

As Fig.|3]and|4]show, the planetary distribution of mass 
and semimajor axis is strongly dependent on the gas ac- 
cretion model considered (Paper I), as well as the type I 
migration rate used. 

When Cmigi = 1 most of the planets migrate to the 
inner disc radius. The mainly difference between the dis- 
tribution found with our model and the one obtained with 



IL04's gas accretion rate (Fig. 3(a) md 4(a) respectively), is 



the population of planets with masses between 1 and IOMq , 
which is larger on the first one. 

In Fig 3(b) and 4(b) {Cmigi = 0.1), the type I migration 
is slower and there are less planets who reach the inner edge 
of the disc. We can see a major population of planets with 
1 — IOM0 and a ^ law. The population of giant planets are 
also increased, but this is larger in Fig |4(b)| because with 
IL04's gas accretion rate the run away of gas is reached 
sooner, and also the Mgap- This allows the planets to start 
the type II migration, which is slower than the type I. 

Type I migration is the slowest when Cmigi ~ 0.01, 
and as seen in Fig 3(c)|4(c)|3(d) and 4(d) the distribution 
found when Cmigi ~ 0.01 is very similar to those found with 
type II migration only. When Cmigi ~ 0.01 the population 
of planets with masses less than 10 is larger than in the 
other cases and there are more giant planets, specially when 
a ^ lau. The differences between 3(c) and 4(c) are similar 
to those found when Cmigi ~ 0.1, which is in agreement 
with the results found in Paper I. 

The "planetary desert" , which is a region with a deficit 
of planets, was found to be between 100 — IOOOM0 (Paper 
I) with our model. Here we found that the effect of the plan- 
etary migration and specially type I migration, is to enlarge 
the desert, allowing to the planets with masses between 10 
and IOOM9 to reach the inner limit of the disc and empty 
the area with semimajor axis between ~ 0.2 and ~ 3au. 
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Figure 3. Distribution of masses and semimajor axis obtained 
with our model. Figur c|3(a) [ shows the data obtained with Ci = 1, 
in [3(b)] Ci = 0.1, m \3(cj\ Ci = .01 and the effect of type I 
migration is not considered in |3(d)| We find a larger population 
of terrestrial planets and a lower population of gas giant beyond 
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Figure 4. Similar figures as [S] but tlie gas accretion rates of 
IL04 are considered. A comparison witli Fig. |3] shows a smaller 
population of terrestrial planets and a larger of giant planets with 
Mt ^ lOOMg), especially when the migration slows down |4(b)] 
[4(^ and [4(d)] 



Table 1. Percentage of final planets tha t are habitable (0.3 ^ 
Mt ^ 7M0 according to llkoma &: Gendal l[2006)'s criteria) found 
with our model (Model I), compared with the results obtained 
with IL04's gas accretion rate (Model II). Different delaying type 
I migration constants are considered. 



Scenario 




Habitable Planets 






{0.3 ^Mt s:7A%) 


Model I, Cmigl 


= 1 


87.02 (%) 


Model II, Cmigl 


= 1 


82.2 (%) 


Model I, C^igi 


= 0.1 


90.5 (%) 


Model II, Cmigl 


= 0.1 


86.9 (%) 


Model I, Cmigl 


= 0.01 


91.57 (%) 


Model II, Cmigl 


= 0.01 


87.18 (%) 



Another important result is that we found a major pop- 
ulation of habitable planets with our model than the one 
found with IL04's gas accretion rates. In order to compare 
the different population of terrestrial planets found with dif- 
ferent gas accretion rates, we adopt a criteria where a hab- 
itable pla net is a planet with m asses between ~ 0.3 and 
« 7M® (|lkoma fc Gendal |2006|). The u pper limit to the 
mass was obtained by Ikoma fc Gendal (|2006| ) with a dif- 



ferent time-scale for gas accretion, with our time-scale this 
value is probably larger, but we consider 7M^ as a nominal 
value for comparison. 

The Table [1] shows th e percentage of final pla nets that 
are habitable according to llkoma fc Gendal l|2006l )'s defini- 
tion. In the table, model I is the one considering our gas 
accretion rate, and model II consider IL04's gas accretion 
rate. 

With our model the run away gas accretion process is 
reached by the cores at a larger mass (Paper I) so there are 
terrestrial planets that never reach this "crossover mass" 
and for this reason we found more planets with masses less 
than 7M0. So a slower gas accretion rate leads to larger 
population of habitable planets. 



4 CONCLUSIONS 

In our previous work, we had developed a very simple model 
for computing planetary formation, which had allowed us 
to show the strong dependence between the gas accretion 
model considered and the mass distribution of extrasolar 
planets. 

In order to get a better understanding of the processes 
involved in planetary formation, we improve our previous 
model, including significant effects that had not been taken 
into account in Paper I while maintaining the computational 
speed. 

We have presented the results of the simulations per- 
formed with our improved model, considering the formation 
of several planets per disc and taking into account the colli- 
sion among them as a source of potential growth. The forma- 
tion of several cores simultaneously in the disc has a strong 
infiuence on the dynamics of the planetesimal disc. The evo- 
lution due to this effect and due to the gas drag effect was 
also considered in this work. 
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Finally, we have analysed the interaction between the 
protoplanets and the disc, which leads to a planetary migrar 
tion (type I and II) . The migration of a giant planet towards 
the central star could perturbate the cores placed in its path, 
causing the ejection or the accretion of the core by the giant 
planet. The another possibility is the survival of the core to 
the passage of the giant planet. We have performed numeri- 
cal simulations taking into account that in a close encounter 
a terrestrial planet has the 71 percent chance of surviving 
the passage of a giant, but these simulations have not pro- 
duced significant changes in the results due to the fact that 
in a real system several protogiant planets cross the region 
of terrestrial planets. 

The final number of planets formed per disc is a direct 
result of the planetary system evolution. We showed that 
the distribution is strongly dependent on the type I migra- 
tion rate considered, noting an absence of planetary systems 
with an intermediate number of planets when the fastest mi- 
gration rate is considered and seeing that it becomes more 
uniform when the migration slowed down. 

The final number of planets per planetary system also 
depends on the initial mass of the disc, those discs with low 
mass will form planetary systems with a large rmmbcr of 
planets and a higher mass disc is necessary to form plane- 
tary systems with a few ones, but the range of disc masses 
corresponding to the final number of planets depends on the 
type I migration rate considered. For a rapid migration rate 
low mass discs arc able to form planetary systems with few 
planets, but when the migration slows down only high mass 
discs could form this kind of planetary systems. 

However we note that many of the systems with a large 
number of planets will be unstable and undergo further evo- 
lution since the orbital spacing is small, and eccentricities 
will increase after the gas is gone. The final number of plan- 
ets will be substantially smaller in some systems, but we do 
not consider this fact in this work. 

We also have analysed the kind of planets formed per 
disc, and found that planetary systems with a large number 
of planets have only terrestrial ones, but when we looked at 
those planetary systems with few planets, the kind of planets 
formed varies according to the type I migration rate consid- 
ered. When we used the fastest rate, most of the planets are 
terrestrial planets, but when the rate is slower the amount 
of terrestrial planets decreases and most are now giants. In 
conclusion we found that the gas giant formation is favored 
at slower type I migration rate. On the other hand the giant 
planets with few gas reminds little in all the simulations, 
this population is not affected by the rate of migration. 

When we analyse the mass and somimajor axis distribu- 
tion we consider different gas accretion rates, and the results 
show that this distribution is strongly dependent on the gas 
accretion model and also on the rate of migration, mainly 
due to type I migration effects. The boundaries of the plan- 
etary desert are enlarged duo to the rapid migration of the 
embryos that did not roach the necessary mass to open up 
a gap. But we note that the planetary desert would be a bit 
different, since the actual observational sample shows that 
the number of observed planets at distances smaller than 
0.1 AU from the central star is not as large as that found 
through numerical simulations, which could mean that most 
planets do not stop their migration at the inner edge of the 
disc. Finally we found that a lower gas accretion rate leads 



to a larger population of habitable planets, and lower pop- 
ulation of gas giants beyond the ice condensation line. 
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